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6  Abstract 

A  mathematical  non-melanoma  skin  cancer  incidence  model  is  derived  which  includes 
the  effects  of  differential  exposure  and  differential  sensitivity  in  the  white 
population  of  a  given  geographic  region.  These  effects,  and  their  variations  with 
geographic  region,  are  shown  to  be  fundamental  to  an  understanding  of  the  behavior 
of  skin  cancer  incidence  as  it  exists  today  and  as  it  would  be  predicted  to  increase 
in  the  event  of  a  decrease  in  stratospheric  ozone.  The  model  is  applied  to  available 
New  York  City  data  to  determine  the  representation  in  the  population  of  nine  postu¬ 
lated  exposure-sensitivity  groups.  The  effects  of  changes  in  exposure  and  sensitiv¬ 
ity  parameters  on  incidence  are  illustrated  by  comparing  New  York  City  with  a  hypo¬ 
thetical  rural  region  at  the  same  latitude.  The  effect  of  not  eliminating  non-solar- 
related  skin  cancers  in  skin  cancer  incidence  data  on  the  prediction  of  increased 
incidence  due  to  stratospheric  ozone  reduction  is  found  to  introduce  an  error  which 
may  be  positive  or  negative.  Programs  which  should  be  pursued  in  order  to  effec¬ 
tively  utilize  the  mathematical  model  are  recommended.  Areas  requiring  urgent  atten¬ 
tion  include  UV-^  instrumentation,  skin  cancer  pathology,  skin  cancer  data  collection 
and  classification,  and  population  exposure  surveys. 
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ABSTRACT 

A  mathematical  non-melanoma  skin  cancer  incidence  model  is 
derived  which  includes  the  effects  of  differential  exposure  and 
differential  sensitivity  in  the  white  population  of  a  given  geo¬ 
graphic  region.  These  effects,  and  their  variations  with  geo¬ 
graphic  region,  are  shown  to  be  fundamental  to  an  understanding 
of  the  behavior  of  skin  cancer  incidence  as  it  exists  today  and 
as  it  would  be  predicted  to  increase  in  the  event  of  a  decrease 
in  stratospheric  ozone.  The  model  is  applied  to  available  New 
York  City  data  to  determine  the  representation  in  the  population 
of  nine  postulated  exposure-sensitivity  groups.  The  effects  of 
changes  in  exposure  and  sensitivity  parameters  on  incidence  are 
illustrated  by  comparing  New  York  City  with  a  hypothetical  rural 
region  at  the  same  latitude.  The  effect  of  not  eliminating  non¬ 
solar-related  skin  cancers  in  skin  cancer  incidence  data  on  the 
prediction  of  increased  incidence  due  to  stratospheric  ozone  re¬ 
duction  is  found  to  introduce  an  error  which  may  be  positive  or 
negative.  Programs  which  should  be  pursued  in  order  to  effec¬ 
tively  utilize  the  mathematical  model  are  recommended.  Areas 
requiring  urgent  attention  include  UV-B  instrumentation,  skin 
cancer  pathology,  skin  cancer  data  collection  and  classification, 
and  population  exposure  surveys. 
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I. 


INTRODUCTION 


Non-melanoma  skin  cancer  Incidence  in  a  given  geographic 
region  is  today  recognized  as  being  strongly  related  to  the  in¬ 
tensity  of  the  solar  ultraviolet  radiation  received,  which  is 
highest  near  the  equator.  Thus,  if  skin  cancer  incidence  is 
plotted  as  a  function  of  latitude,  it  would  be  expected  that 
there  would  be  a  well-defined  latitude  gradient.  In  Fig.  1  is 
reproduced  a  plot  of  non-melanoma  skin  cancer  incidence  for  male 
white  populations  between  30°  and  60°  north  latitude  (Ref.  1). 
While  there  is  a  discernible  latitude  gradient  in  this  plot, 
there  is  a  very  large  scatter  between  40°  and  55°  which  can  be 
principally  attributed  to  the  combined  effects  of  poor  data, 
differences  in  climate,  differences  in  ozone  thickness,  differ¬ 
ences  in  sun  exposure  habits  of  the  population  (which  are  related 
to  climate),  and  differences  in  population  sensitivity  to  non¬ 
melanoma  skin  cancer.  The  complex  ways  in  which  the  latter  two 
effects  influence  skin  cancer  incidence  are  explored  In  this 
paper.  A  simplified  exposure-sensitivity  model  is  derived  and 
applied  to  available  data  for  New  York  City.  A  comparison  is 
then  made  with  a  hypothetical  rural  region  of  the  same  latitude, 
in  which  exposure  and  sensitivity  parameters  can  be  expected  to 
differ  widely. 

Nearly  all  of  the  non-melanoma  skin  cancer  incidence  models 
that  have  been  proposed  for  predicting  the  increase  in  incidence 
resulting  from  an  increase  in  solar  ultraviolet  radiation  due  to 


a  postulated  decrease  In  stratospheric  ozone  are  aggregate  models, 
i.e..  It  Is  implicitly  assumed  that  all  members  of  a  given  white 
population  are  exposed  to  the  same  amount  of  UV-B  solar  radiation 
(e.g.,  some  fraction  of  the  mean  UV-B  radiation  incident  on  the 
ground)  and  have  the  same  sensitivity  to  skin  cancer.  In  actu¬ 
ality,  of  course,  there  are  very  wide  differences  in  both  the  sun 
exposure  habits  and  sensitivity  of  individuals.  Consequently,  it 
can  be  expected  that  individuals  who  seldom  expose  themselves  to 
the  sun  and  who  have  low  sensitivity  to  skin  cancer  would  be  less 
affected  by  a  postulated  decrease  in  the  amount  of  stratospheric 
ozone  than  those  individuals  who  have  high  sensitivity  and  high 
sun  exposure.  However,  if  this  high-exposure,  high-sensitivity 
group  constitutes  but  a  small  fraction  of  the  overall  population, 
its  contribution  to  the  value  of  skin  cancer  incidence  for  the 
whole  population  could  be  small.  It  also  follows  that  any  addi¬ 
tional  cases  to  be  expected  from  this  high  risk  group  as  a  re¬ 
sult  of  a  postualted  decrease  in  stratospheric  ozone  could  con¬ 
stitute  but  a  small  fraction  of  the  total  additional  cases.  In 
the  model  derived  below,  skin  cancer  incidence  is  derived  as  a 
function  of  the  probability  distributions  of  exposure  and  sensi¬ 
tivity  in  the  population.  To  adequately  describe  these  distri¬ 
butions,  it  is  necessary  to  assume  ranges  of  relative  exposure 
and  sensitivity  values  and  the  probability  functions  within  these 
assumed  ranges.  The  model  indicates  that  the  skin  cancer  inci¬ 
dence  following  any  calculated  change  in  the  ozone  column,  e.g., 
based  on  2-D  modelling  of  the  effects  of  aircraft  effluents, 
will  depend  significantly  on  the  differential  sensitivity  and 
exposure  characteristics  of  the  population  residing  in  a  given 
geographic  region,  as  well  as  on  the  magnitude  of  the  change  in 
the  ozone  column. 
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II.  MODEL  FOR  SKIN  CANCER  INCIDENCE 


Non-melanoma  skin  cancer  incidence  increases  very  rapidly 
with  age.  When  plotted  on  semi-log  paper,  age-specific  non¬ 
melanoma  skin  cancer  incidence  typically  approximates  a  straight 
line,  as  indicated  by  the  diverse  smoothed  data  in  Fig.  2  for 
incidence  in  Finland,  Sweden,  Japan  (Ref.  2)  and  mortality  in 
the  U.S.  (Ref.  3).  The  slopes  of  the  smoothed  lines  are  almost 
equal;  their  values  indicate  an  incidence  doubling  time  of  6  to 
7  years.  For  a  population  group  of  age  A,  exposure  E,  and  sen¬ 
sitivity  S,  the  incidence  I'  can  therefore  be  represented  approx¬ 
imately  by 


I'  (A,  S,  E)  =  k  f  (S,  E)  eXA,  (1) 

where  k  and  A  are  constants. 

The  term  "exposure,"  here  designated  by  E,  represents  the 
received  biologically  weighted  UV-B  (sometimes  referred  to  as 
DUV  or  damaging  UV)  flux  integrated  over  exposure  time.  It  is 
assumed  here  that  the  biological  weighting  function  (usually 
assumed  to  be  the  erythema  action  spectrum)  has  the  same  shape 
for  all  individuals.  However,  Caucasians  are  known  to  differ 
in  their  tendency  to  develop  solar-related  skin  cancer  because 
of  differences  in  the  UV  reflectance  and  absorption  properties 
of  the  human  skin  and  differences  in  DNA  repair  mechanisms;  this 
complex  biological  difference  is  covered  by  the  sensitivity 
factor  S. 

The  function  f  (S,  E)  is  assumed  to  follow  the  power  law 
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FIGURE  2.  Non-melanoma  skin  cancer  incidence  and 
mortality  vs  age  (heavily  smoothed). 
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(2) 


f  (S,  E)  =  (SE)n  =  Dn  . 

The  product  SE  in  Eq.  (2)  is  a  measure  of  the  carcinogenically 
effective  received  ultraviolet  dose  D  for  a  population  of  expo¬ 
sure  E  and  sensitivity  S,  and  n  is  the  biological  amplification 
factor,  i.e.. 


dl/I 
n  “  dD/D 


(3) 


n  is  apparently  independent  of  age  A,  since  the  distance  between 
any  two  lines  for  two  different  received  dose  levels  in  Fig.  2 
is  essentially  a  constant.  The  higher  the  value  of  n,  of  course, 
the  higher  is  the  percentage  increase  in  skin  cancer  incidence 
for  a  specified  percentage  increase  in  UV  dose. 


If  p ( A) ,  p (E) ,  and  p(S)  represent  the  probability  density 
functions  for  population  age,  exposure,  and  sensitivity,  respec¬ 
tively,  incidence  I  for  the  whole  population  of  a  given  region 
is  given  by 


I 


n  eXA  p(S)ps(E)p(A)dSdEdA 
(S)ps(E)dSdE, 


(4) 

(5) 


where  p„(E)  is  the  conditional  probability  of  exposure  E  given 
a  sensitivity  S,  and  k'  is  a  constant. 

In  order  to  reduce  the  triple  integral  (4)  to  a  double  in¬ 
tegral  (5)  by  Integrating  over  age,  it  was  implicitly  assumed 
that  p(S)  and  Pg(E)  are  independent  of  age.  For  a  stable  popu¬ 
lation  composition  such  as  that  to  be  found  in  Sweden,  p(S)  in¬ 
deed  should  be  close  to  being  independent  of  age.  However, 
p  (E)  for  adults  could  be  expected  to  be  dependent  on  age  as 
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well  as  sensitivity.  Older  retired  people  who  are  healthy  prob¬ 
ably  spend  more  time  outdoors  than  they  did  when  working,  while 
the  infirm  probably  spend  less  time  outdoors.  While  these  ef¬ 
fects  tend  to  offset  each  other,  the  extent  to  which  they  do  not 
represents  a  deficiency  in  this  representation. 

The  continuous  probability  distributions  in  Eq.  (5)  are  un¬ 
known  quantities  at  the  present  time.  It  is  easier  to  work  with 
a  simplified  discrete  form  of  Eq.  (5),  i.e., 

i=a  j=b 

1  ■  k'  2  Z  (Ei  Vn  Py  >  (6) 

i=l  j=l 

where  each  population  subgroup  of  exposure  E.^  and  sensitivity 
has  a  probability  p^  of  being  represented  in  the  total  popula¬ 
tion.  The  latter  is  divided  into  "a"  exposure  subgroups  and  "b" 
sensitivity  subgroups,  and 


i=a  j  =b 

Z  • 

i=l  j-1 


III.  ANALYSIS  OF  SKIN  CANCER  INCIDENCE  IN  NEW  YORK  CITY 


Very  little  quantitative  information  is  available  in  the 
literature  comparing  incidences  of  population  subgroups  of  dif¬ 
ferent  exposures  and  sensitivities  in  a  given  geographical  re¬ 
gion.  The  only  quantitative  information  of  this  kind  known  to 
the  author  is  for  a  basal  cell  epithelioma  control  study  for 
New  York  City  (Ref.  4).  The  data  in  this  study  fortuitously 
provides  values  for  some  of  the  exposure  and  sensitivity  param¬ 
eters  required  in  the  analysis  below.  It  should  be  noted  that 
while  the  bulk  of  non-melanoma  skin  cancers  are  basal  cell  epi¬ 
theliomas,  squamous  cell  carcinomas  constitute  a  sizeable  frac¬ 
tion  of  all  skin  cancers.  Moreover,  not  all  basal  cell  epithe¬ 
liomas  are  induced  by  solar  ultraviolet  radiation  (Ref.  1). 

Gellin  et  al.  (Ref.  4)  divided  the  population  into  three 
exposure  subgroups,  according  to  the  estimated  daily  hours  of 
outdoor  activity.  The  low-exposure  subgroup  spent  0  to  2  hours 
daily  outdoors,  the  moderate-exposure  group  3  to  5  hours,  and 
the  high-exposure  group  6  or  more  hours.  The  percentage  distri¬ 
bution  of  these  subgroup  populations  in  the  BCE  (basal  cell 
epithelioma)  and  control  groups  are  given  in  Table  1  (Ref.  4). 
The  data  set  used  consisted  of  759  cases  of  basal  cell  epithe¬ 
lioma  and  757  controls,  and  was  restricted  to  Caucasians. 

TABLE  1.  PERCENTAGE  DISTRIBUTION  BY  ESTIMATED  DAILY 
HOURS  OF  OUTDOOR  ACTIVITY 


EXPOSURE  (hrs ) 

BCE 

CONTROLS 

0-2 

24.7 

65.0 

3-5 

46.4 

25.1 

6  + 

28.9 

9.9 

9 


Information  on  a  highly  sensitive  subgroup  characterized 
by  light  hair,  fair  complexion,  and  light  eyes  was  also  reported, 
independent  of  exposure,  and  for  that  portion  of  this  particular 
subgroup  that  estimated  an  outdoor  exposure  of  6  or  more  hours 
daily.  The  percentage  distributions  are  given  in  Table  2  for 
the  BCE  and  control  groups  (Ref.  4) . 

TABLE  2.  PERCENTAGE  DISTRIBUTION  FOR  HIGHLY  SENSITIVE  SUBGROUPS 

SUBGROUPS  BCE  CONTROLS 

Highly  Sensitive  15.3  4.2 

Highly  Sensitive  and  5.3  0.6 

6+  Hrs .  Outdoor  Activity 


An  inspection  of  Tables  1  and  2  reveals  the  significance  of 
outdoor  exposure  on  skin  cancer  incidence,  particularly  for  those 
who  are  highly  sensitive. 


The  New  York  City  population  is  assumed  to  be  divided  into 
the  three  exposure  groups:  E1,  E2,  and  E^,  as  indicated  in 
Table  1,  and  three  sensitivity  groups:  S^,  S2,  and  S^.  The 
subscript  1  denotes  a  low,  2  a  moderate,  and  3  a  high  value  of 
exposure  or  sensitivity.  Consequently,  there  will  be  3  x  3  or 
nine  sensitivity-exposure  groups,  each  characterized  by  a  dose 
(E^  Sj)n  as  indicated  in  the  matrix  of  Fig.  3  for  n  *  1. 


E1 

EXPOSURE 

E2 

E3 

S1 

E,S, 

E2S1 

E3S1 

S2 

E1S2 

E2S2 

E2S2 

S3 

E1S3 

E2S3 

E3S3 

FIGURE  3.  DOSE  MATRIX  (n*l) 
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For  each  of  the  population  subgroups  in  Fig.  3  characterized 
by  dose  Sj ,  there  is  associated  a  probability  p^j ,  as  indicated 
by  the  probability  distribution  matrix  of  Fig.  4.  The  probability 
p^j  may  be  interpreted  as  either  the  probability  of  finding  a  ran¬ 
domly  chosen  member  of  the  whole  population  who  has  exposure 
and  sensitivity  Sj ,  or  the  fraction  of  the  whole  population  who 
have  exposure  E.^  and  sensitivity  Sj  . 


OO 

Z 

Ui 

go 


E1 

EXPOSURE 

E2 

E3 

si 

pn 

P21 

P31 

S2 

p12 

P22 

P32 

S3 

PI  3 

P23 

P33 

FIGURE  4.  PROBABILITY  DISTRIBUTION  MATRIX 


From  Table  1  and  Fig.  it  follows  that 

P11  +  P12  +  p13  =  0,65  (7) 

P21  +  P22  +  p23  =  0,25  (8) 

P31  +  P32  +  p33  =  0,10  (9) 

If  f^,  f2>  and  f^  represent  the  fractions  of  the  New  York 

City  population  with  sensitivities  S^,  S2,  and  respectively, 
then  from  Fig.  1| 


P11  +  P21  +  P31  =  fl 
P12  +  P22  +  P32  *  f2 
P13  +  p23  +  P33  f3  1_fl”f2 
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since  the  f^  values  are  subject  to  the  constraint 


i-3 


L 

4  si 


1 


(13) 


Eq.  (12)  is  not  an  independent  equation  since  it  can  be  derived 
from  Eqs .  (7),  (8),  (9),  (10),  and  (11). 

Prom  Eq.  (6),  Table  1,  Fig.  1,  and  Fig.  3  it  follows  that 
the  ratio  of  the  contribution  to  BCE  incidence  from  the  moderate 
exposure  group  (column  2  of  each  matrix)  to  the  contribution 
from  the  low  exposure  group  (column  1  of  each  matrix)  is  given 
by 


E2  S1  P11  *  E2  S2  p22  +  E2  S3  P23  46.4 

E1  S1  P11  +  E1  S2  p12  +  E1  S3  P13  2k’7 


=  1.88 


(14) 


and  that  the  ratio  of  the  contribution  to  BCE  incidence  from  the 
high  exposure  group  (column  3  of  each  matrix)  to  the  contribution 
from  the  low  exposure  group  is  given  by 


E3  S1  P31  *  E3  S2  p32  *  E3  S3  P33  =  28.9 
E1  S1  P11  +  E1  ^2  p12  +  E1  S3  P13  2 21,7 

=  1.17 


(15) 


From  Table  2  and  Fig.  4  it  can  be  seen  that  the  ratio  of 
the  percentage  of  New  York  City  people  who  are  in  the  highly  sen¬ 
sitive  group  to  the  percentage  that  are  highly  sensitive  and 
spend  6  or  more  hours  daily  outdoors  is  given  by 
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(16) 


P13  *  P23  +  p33  _  4.2  _ 
p33  -  OT^  “ 

or 

^23  ~  6,0  “0  •  (17) 

From  Table  2,  Eq.  (6),  Fig.  3,  and  Fig.  4  it  can  be  seen 
that  the  ratio  of  the  percentage  of  BCE  cases  in  the  highly 
sensitive  group  to  the  percentage  of  BCE  cases  in  the  highly 
sensitive  group  who  also  spend  6  or  more  hours  daily  outdoors 
is  given  by 


E1  S3  P13  *  E2  S3_P23  *  E3  S3  p33  =  16-3 


E3  S3  p33 


5.3 


(18) 


or 

El  p13  +  E2  p23  “  1*887  E3  P33  =  0  •  d9) 

If  values  for  f^  f2,  E1,  E2>  E^,  S1,  S2,  and  are 
assigned,  it  is  possible  to  solve  the  nine  independent  linear 
equations  7,  8,  9,  10,  11,  14,  15,  17,  and  19  for  the  nine  p^ 
values.  Since  only  relative  values  of  exposure  and  sensitivity 
are  necessary  inputs  to  the  solution,  it  Is  convenient  to  assign 
a  value  of  unity  for  the  low  exposure  and  sensitivity  groups, 
i.e . , 

E1  =  1  (20) 

Sx  -  1  (21) 

The  remaining  two  values  of  group  exposure  and  sensitivity, 
by  definition,  must  satisfy  the  inequalities 

1  <  E2  <  E3  (22) 

1  <  S2  <  S3  .  (23) 


13 


The  values  of  and  determine  the  range  of  exposure  and 
sensitivity.  Unfortunately,  there  is  no  quantitative  information 
available  on  sensitivity  variations  within  white  populations. 
There  is  some  guidance  for  the  selection  of  E  values  in  the 
numerical  values  of  estimated  daily  hours  of  exposure  in  Table 
1.  However,  there  are  the  following  difficulties  in  selecting 
appropriate  E  values  for  group  exposure: 

1.  The  low  exposure  group  spends  0  to  2  hours  outdoors 
daily,  and  it  is  tempting  to  take  a  mean  value  of 
one  hour  as  representative  of  this  group,  but  the 
distribution  within  this  2-hour  period  is  unknown; 
the  time  mean  could  lie  significantly  below  or  above 
the  1-hour  value. 

2.  Time  spent  outdoors  is  not  a  very  good  measure  of  the 
received  solar  ultraviolet  dose,  but  it  is  the  only 
data  available.  Exposures  during  the  hours  10:00  AM  to 
2:00  PM,  of  course,  result  in  much  higher  UV  doses  than 
early  morning  or  late  afternoon  exposures  of  equal  dura¬ 
tion.  Summer  exposures  also  result  in  higher  UV  doses 
than  exposures  during  the  other  seasons  for  the  same 
time  period. 

There  exists  an  acceptability  constraint  on  the  set  of  the 
input  values  E2,  E^,  S2,  S^,  f^,  and  f2  that  are  selected:  the 
solution  of  the  simultaneous  equations  must  result  in  a  set  of 
nine  positive  p^j  values,  i.e.,  a  negative  value  of  probability 
has  no  physical  meaning. 

Parameter  values  are  listed  in  Table  3  for  22  cases  which 
were  selected  for  solution  by  a  computer.  Of  these,  seven  cases 
were  found  to  result  in  a  solution  with  all  nine  p^  values 
positive  (/);  the  remaining  15  cases  resulted  in  one  or  more 
negative  values  (X)  of  p^  and  hence  are  unacceptable. 
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w" 

? 

1. 

Case 

n 

/ 

TABLE  3. 

E2 

ACCEPTABILITY  OF 

E3  S2  S3 

CASES 

fl 

1 

INVESTIGATED  1 

f g  Acceptability  I 

A 

1 

4 

8 

2 

4 

0.2 

0.6 

/ 

t- 

B 

1 

4 

8 

2 

4 

0.1 

0.8 

X 

[■ 

C 

1 

4 

8 

2 

4 

0.3 

0.4 

X 

P-. 

0 

1 

4 

8 

2 

3 

0.2 

0.6 

/ 

» 

■ 

E 

1 

4 

8 

2 

3 

0.1 

0.8 

X 

L 

F 

1 

4 

8 

2 

3 

0.3 

0.4 

X 

\ 

G 

1 

4 

8 

3 

9 

0.2 

0.6 

/ 

i 

1  H 

i 

1 

4 

8 

3 

9 

0.1 

0.8 

/ 

I 

1 

4 

8 

3 

9 

0.3 

0.4 

X 

J 

1 

4 

8 

3 

9 

0.5 

0.3 

X 

, 

K 

1 

4 

8 

3 

9 

0.2 

0.3 

X 

1  L 

1 

4 

6 

2 

4 

0.2 

0.6 

X 

M 

1 

4 

6 

2 

4 

0.1 

0.8 

X 

i  N 

1 

4 

6 

2 

4 

0.3 

0.4 

/ 

0 

1 

4 

6 

2 

4 

0.4 

0.2 

/ 

f 

P 

1 

4 

6 

3 

6 

0.2 

0.6 

X 

| 

Q 

1 

4 

6 

3 

6 

0.1 

0.8 

X 

1  R 

I 

1 

4 

6 

3 

6 

0.3 

0.4 

X 

t 

s 

2 

4 

8 

2 

3 

0.2 

0.6 

X 

T 

2 

4 

8 

2 

3 

0.1 

0.8 

X 

U 

2 

4 

8 

2 

3 

0.3 

0.4 

X 

V 

2 

4 

8 

2 

3 

0.15 

0.7 

/ 

15 


/ 


The  p. .  matrices  for  the  seven  acceptable  cases  listed  in 
Table  3  are  given  in  Fig.  5.  All  cases  except  the  last  four 
assumed  a  linear  dose-response  relationship,  i.e.,  n  =  1.  All 
cases  except  unacceptable  cases  J  and  K  assumed  a  symmetric 
distribution  of  sensitivity  in  the  population,  i.e.,  f^  *  f^ 
or  f2  *  1  -  2  f1»  For  each  of  five  sets  of  exposure  and  sen¬ 
sitivity  values,  the  case  f^  =  0.2  was  assumed,  and  additional 
cases  below  and  above  this  value  were  investigated  to  determine 
how  sensitive  the  solution  was  to  this  parameter.  If  all  the 
parameters  except  f1  are  held  constant,  each  p  is  a  linear 
function  of  f^,  as  illustrated  in  Fig.  6,  which  encompasses  cases 
A,  B,  and  C.  Note  that  in  order  for  all  p^  values  to  be  posi¬ 
tive,  the  particular  conditions  cited  in  Fig.  5  require  that 
0.134  <  ^  <  0.240. 

The  seven  acceptable  solutions  in  Table  3  can  not  all  be 
correct;  Indeed,  it  cannot  be  expected  that  any  one  of  these  so¬ 
lutions  could  be  correct  for  the  following  principal  reasons: 

1.  The  set  of  cases  investigated  is  only  a  small  fraction 
of  the  possible  number  of  sets  that  could  be  obtained 
by  taking  different  combinations  of  the  assumed  param¬ 
eter  values  in  Table  3  and  by  considering  additional 
parameter  values. 

2.  The  model  is  extremely  crude,  e.g.,  the  assumption  of 
three  exposure  and  three  sensitivity  groups  is  a  gross 
simplification  of  the  physical  situation  involving  con¬ 
tinuous  distributions  in  both  exposure  and  sensitivity. 

3.  The  input  incidence  data  for  New  York  City  is  crude, 
since  it  is  based  on  a  total  of  only  759  BCE  cases  In 
the  period  January  1955  to  March  1959  (Ref.  4). 

4.  New  York  City  has  a  polyglot  population  with  a  large 
fraction  of  immigrants  who  experienced  a  different  UV 
environment  in  their  native  lands.  Non-melanoma  skin 
cancer  incidence  is  a  function  of  lifetime  UV  dose  (Ref. 5). 
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.1223 

.0097 
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.1017 

.0571 

PROBABILITY  MATRICES  FOR 
SEVEN  ACCEPTABLE  CASES 
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5.  The  input  incidence  data  used  for  New  York  City  was 
for  the  combined  sexes.  Exposure  behavior  and  cloth¬ 
ing,  and  consequently  incidence,  vary  significantly 
with  sex. 

6.  Perhaps  as  many  as  one-half  of  basal  cell  carcinomas 
are  not  induced  by  solar  ultraviolet  radiation  in 
northern  U.S.  cities  (Ref.  1). 

7.  Squamous  cell  carcinomas  were  not  included. 

8.  Differences  in  exposure  with  age  were  not  considered. 

Any  future  models  including  age  as  a  parameter  will 
require  three-dimensional  matrices. 

It  was  deemed  not  worthwhile  to  investigate  additional 
cases  in  response  to  the  first  reason  cited  above  because  of 
the  limitations  imposed  by  the  other  reasons.  It  is  possible, 
however,  to  devise  criteria  which  enable  comparisons  to  be  made 
between  the  seven  acceptable  solutions  in  order  to  select  the 
solution  that  best  fits  the  incidence  data. 

In  Table  4  the  probability  ratios  are  given  for  groups  of 
equal  sensitivity  and  different  exposure.  Thus,  for  example, 
the  ratio  is  the  ratio  °T  the  fraction  of  the  least  sen¬ 

sitive  population  which  spends  0  to  2  hours  daily  outdoors  to 
the  fraction  which  spends  3  to  5  hours  daily  outdoors.  For  the 
entire  population,  the  ratio  of  the  fraction  spending  0  to  2 
hours  daily  outdoors  to  that  spending  3  to  5  daily  outdoors  is 
0.65/0.25  (Table  1)  or  2.6. 

One  would  expect  that  the  ratio  of  time  spent  outdoors  for 
the  sensitive  group  to  be  greater,  but  not  very  much  greater, 
than  the  ratio  for  the  entire  population.  The  ratio  of  3.2  for 
Case  A  in  Table  4  fulfills  this  expectation.  Similarly,  the 
ratio  P23/P31  T°r  the  Traction  of  the  least  sensitive  group  which 
spends  3  to  5  hours  outdoors  daily  to  the  fraction  spending  more 
than  6  hours  outdoors  daily  would  not  be  expected  to  deviate 
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TABLE  4.  PROBABILITY  RATIOS  FOR  GROUPS  OF  EQUAL  SENSITIVITY 


greatly  from  0.25/0.10  (Table  1)  or  2.5.  The  ratio  of  0.6  for 
Case  A  in  Table  4  is  one-fourth  of  that  for  the  entire  popula¬ 
tion,  which  seems  unlikely  but  is  not  impossible  since  the  least 
sensitive  group  would  be  more  likely  to  seek  high  exposure  than 
the  population  as  a  whole. 

The  second  and  third  columns  in  Table  4  are  the  ratios  for 
the  moderate  sensitivity  group,  and  the  fifth  and  sixth  columns 
are  the  ratios  for  the  high  sensitivity  group.  Thus,  for  a 
given  case,  there  are  six  ratios  to  be  compared,  and  all  six  ra¬ 
tios  would  not  be  expected  to  deviate  greatly  from  a  value  of  2.5. 

As  a  measure  of  the  spread  in  the  values  of  the  ratios  for 
a  given  case,  the  maximum  value  is  divided  by  the  minimum  value 
(last  column  in  Table  4).  A  large  max  ratio/min  ratio  value 
indicates  a  large  excursion  from  expected  behavior.  This  value 
fell  in  the  range  of  15  to  18  for  Cases  A,  G,  H,  and  0,  but  ex¬ 
ceeded  120  for  Cases  D,  N,  and  V.  The  latter  three  cases  repre¬ 
sent  too  large  a  departure  from  expected  behavior  and  are  there¬ 
fore  discarded. 

In  Table  5,  the  probability  ratios  are  given  for  groups  of 
equal  exposure  and  different  sensitivity  weighted  by  the  appro¬ 
priate  f  ratio.  Only  the  four  remaining  Cases  A,  G,  H,  and  0 
are  considered.  If  f^  =  f^,  the  ratio  P11/P12  would  represent 
the  ratio  of  the  fraction  of  the  low  exposure  population  group 
of  low  sensitivity  to  the  fraction  of  the  low  exposure  group  of 
moderate  sensitivity.  Since,  in  the  cases  investigated,  f^  ^  f2, 
the  ratio  P11/P12  is  muitiPlled  by  If  exposure  were  in¬ 

dependent  of  sensitivity,  it  would  be  expected  that  P1]_/P12  x 
fj/fi  would  be  equal  to  unity.  Values  in  Table  5  which  deviate 
greatly  from  unity  are  considered  less  likely  to  be  realistic 
than  values  which  are  close  to  unity. 

By  again  applying  the  max  ratio/min  ratio  criterion  to  the 
four  cases  in  Table  5>  a  minimum  value  of  15  is  obtained  for 
Case  H,  and  a  maximum  value  of  50  for  Case  G.  Thus,  the  H  set 
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of  input  values  n=  1,  E2  =  4,  E^  =  8,  S2  =  3,  =  9,  and 

f1  =  =  0.1  appears  to  best  fit  the  input  New  York  City  data 

if  the  population  is  divided  into  nine  discrete  sensitivity- 
exposure  groups  and  the  linear  dose-response  relationship  is 
assumed  to  be  valid.  The  fact  that  an  acceptable  but  non¬ 
competitive  solution  was  found  for  n  =  2  (Case  V  in  Table  3) 
indicates  that  there  is  a  possiblity  that  a  biological  amplifi¬ 
cation  factor  n  other  than  unity  could  lead  to  a  better  fit  of 
the  data. 


The  fraction  1^ j  of  the  total  skin  cancer  incidence  con¬ 
tributed  by  a  group  with  exposure  and  sensitivity  Sj  is 
given  by 


I 


ij 


E1  SJ_  PU 

i=3  j=3 

ZEh  ^  ^ 

i=l  J=1 


(24) 


In  Fig.  7,  the  elements  of  the  dose  (E^  Sj ) ,  probability 
(p^j),  and  incidence  (1^)  matrices  are  given  for  Case  H.  From 
the  incidence  fraction  matrix  it  is  seen  that  while  only  1.4 
percent  of  the  population  is  in  the  high  exposure-high  sensitivity 
group,  12.1  percent  of  skin  cancer  cases  occur  in  this  group. 

On  the  other  hand,  whereas  5  percent  of  the  population  is  in 
the  low  exposure-low  sensitivity  group,  only  0.6  percent  of  the 
skin  cancer  cases  occur  in  that  group.  These  results  are  con¬ 
sistent  with  the  fact  that  a  high  exposure-high  sensitivity 
resident  runs  a  risk  of  skin  cancer  72  times  higher  than  a  low 
exposure-low  sensitivity  resident,  as  indicated  by  the  dose 
matrix  in  Fig.  7,  i.e.,  12.1/1.4  *  .  006/.05  =  72.  While  a  ma¬ 
jority  of  the  population  (55.8  percent)  lies  in  the  low  exposure- 
moderate  sensitivity  group,  only  19.6  percent  of  the  skin  cancer 
cases  are  found  in  this  group.  The  group  contributing  most  skin 
cancers  is  the  moderate  exposure-moderate  sensitivity  group,  with 
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27.6  percent  compared  to  a  population  representation  of  19-6 
percent.  Both  the  high  exposure-low  sensitivity  group  and  the 
low  exposure-high  sensitivity  group  contribute  approximately 
4  percent  to  skin  cancer  incidence  and  have  approximately  the 
same  4  percent  population  representation. 

As  shown  in  the  incidence  fraction  matrix  of  Fig.  7,  ad¬ 
dition  of  each  element  in  a  row  yields  the  incidence  fraction 
for  a  group  of  given  sensitivity,  and  addition  of  each  element 
in  a  column  yields  the  incidence  fraction  for  a  group  of  given 
exposure.  In  Table  6,  the  population  and  incidence  percentages 
are  compared  for  the  three  assumed  levels  of  exposure  and  sen¬ 
sitivity,  i.e.,  low,  moderate,  and  high.  Residents  with  high 
exposure  or  high  sensitivity  contribute  to  incidence  at  a  rate 
approximately  triple  that  of  their  representation  in  the  popu¬ 
lation.  Residents  with  low  exposure  or  low  sensitivity  contrib¬ 
ute  at  a  rate  less  than  half  their  representation  in  the  popu¬ 
lation.  Residents  of  moderate  exposure  contribute  at  a  rate 
almost  double  their  population  representation,  while  those  of 
moderate  sensitivity  contribute  at  a  rate  approximately  three- 
fourths  of  their  population  representation.  Residents  with  high 
exposure  or  high  sensitivity  contribute  at  a  rate  three  times 
higher  than  their  representation  in  the  population. 

TABLE  6.  POPULATION  AND  INCIDENCE  PERCENTAGES  FOR 

NEW  YORK  CITY 

EXPOSURE  SENSITIVITY 


LEVEL 

POPULATION 

INCIDENCE 

POPULATION 

INCIDENCE 

Low 

65 

24.7 

10 

4.8 

Moderate 

25 

46.4 

80 

60.3 

High 

10 

28.9 

10 

34.9 

25 


I 


The  I^j  values  in  Fig.  7  also  give  information  on  the 
sources  of  new  skin  cancer  cases  in  the  event  of  a  reduction 
in  the  amount  of  stratospheric  ozone.  The  number  of  new  cases 
would  be  found  to  have  an  exposure-sensitivity  distribution  that 
would  parallel  the  I^j  distribution,  i.e.,  each  I^j  would  re¬ 
ceive  the  same  percentage  increase  and  hence  any  two  elements 
of  a  matrix  of  number  of  new  skin  cancer  cases  would  have  the 
same  ratio  as  their  1^  ratio. 
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IV.  RURAL  CASE 


It  is  interesting  to  investigate  the  change  in  the  elements 
of  the  incidence  matrix  in  Fig.  7  if  a  hypothetical  rural  region 
were  to  be  considered  instead  of  New  York  City.  In  a  region 
such  as  Iowa,  which  lies  at  approximately  the  same  latitude  as 
New  York  City  and  where  farming  is  the  predominant  occupation, 
it  might  be  expected  that  a  majority  of  the  population  spends 
more  than  6  hours  outdoors  daily*  instead  of  the  10  percent  for 
New  York  City. 

In  the  hypothetical  rural  example  of  Fig.  8,  the  fraction 
of  the  population  spending  6  or  more  hours  outdoors  daily  is 
assumed  to  be  65  percent  and  the  fraction  spending  0  to  2  hours 
outdoors  daily  is  assumed  to  be  10  percent,  or  the  reverse  of 
the  New  York  City  exposure  distribution.  The  sensitivity  dis¬ 
tribution  is  assumed  to  be  the  same  as  that  for  New  York  City, 
i.e.,  f1  *  0.1  and  f2  -  0.8.  The  moderate  sensitivity-high 
exposure  group  was  assumed  to  constitute  55  percent  of  the  pop¬ 
ulation,  a  value  which  corresponds  to  that  of  the  moderate 
aensitivity-low  exposure  group  in  New  York  City.  In  this  rural 
example,  83  percent  of  the  skin  cancers  are  in  the  high  exposure 
group.  A  comparison  of  the  population  and  incidence  percentages 
for  this  hypothetical  rural  case  is  given  in  Table  7.  Note  in 
Table  7  that  the  high  sensitivity  group  accounts  for  only  21 
percent  of  the  Incidence  whereas  in  Table  6  it  accounted  for 
35  percent.  The  explanation  for  this  surprising  result  is  the 
1 - 

Not  likely  to  be  true  for  the  female  population,  but  con¬ 
sideration  of  the  sex  difference  is  not  required  for  pur¬ 
poses  of  comparison  of  this  hypothetical  rural  case  with 
the  data  set  of  Section  III. 


PROBABILITY  (p^) 


j  =  l 


INCIDENCE  FRACTION  (I..) 

*  J 


FIGURE  8.  DOSE,  PROBABILITY,  AND  INCIDENCE  MATRICES  FOR  A 
RURAL  HYPOTHETICAL  EXAMPLE  WITH  A  NEW  YORK  CITY 
SENSITIVITY  DISTRIBUTION 


28 


assumed  order  of  magnitude  increase  in  the  high  exposure-moderate 
sensitivity  population  fraction  where  the  majority  (64  percent) 
of  rural  skin  cancers  are  found. 


TABLE  7.  POPULATION  AND  INCIDENCE  PERCENTAGES  FOR 
A  HYPOTHETICAL  RURAL  REGION  WITH  A 
NEW  YORK  CITY  SENSITIVITY  DISTRIBUTION 


EXPOSURE  SENSITIVITY 


LEVEL 

POPULATION 

INCIDENCE 

POPULATION 

INCIDENCE 

Low 

10 

2.5 

10 

2.8 

Moderate 

25 

14.2 

80 

76.3 

High 

65 

83.3 

10 

20.9 

The  ratio  of  skin  cancer  incidence  of  this  rural  example 
to  the  incidence  for  New  York  City  is  equal  to  the  ratio  of  the 
respective  denominators  of  Eq.  (24),  which  is  equal  to  20.64/ 
8.5194  or  2.42. 

In  Fig.  9  the  effect  of  a  sensitivity  distribution  change 
has  been  investigated  in  which  the  fraction  of  the  high  sensi¬ 
tivity  group  assumed  in  Figs.  7  and  8  is  doubled  from  0.1  to 
0.2.  This  increase  is  assumed  to  be  associated  with  a  corre¬ 
sponding  decrease  in  the  moderate  sensitivity  group  and  no  change 
in  the  low  sensitivity  group,  i.e.,  f^  *  0.1,  f2  =  0.7,  and 
f^  =  0.20.  A  comparison  of  the  population  and  Incidence  per¬ 
centages  for  this  rural  case  is  given  in  Table  8.  The  percent¬ 
ages  by  level  of  exposure  remained  essentially  unchanged  from 
the  percentages  in  Table  7.  The  fraction  of  the  Incidence  con¬ 
tributed  by  the  high  sensitivity  group  is  almost  double  the  pop¬ 
ulation  fraction  of  20  percent.  However,  because  70  percent  of 
the  population  falls  in  the  moderate  sensitivity  group,  the  total 
incidence  for  the  case  in  Fig.  9  increased  by  a  factor  of  only 
23.7/20.64,  or  1.15  over  the  case  in  Fig.  8. 
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FIGURE  9.  DOSE,  PROBABILITY,  AND  INCIDENCE  MATRICES  FOR  A 

RURAL  HYPOTHETICAL  EXAMPLE 
WITH  f,  =  0.1  and  f«  =  0.2 


TABLE  8.  POPULATION  AND  INCIDENCE  PERCENTAGES  FOR 
A  HYPOTHETICAL  RURAL  REGION  WITH 
=  0.1  and  =  0.2 


EXPOSURE  SENSITIVITY 


LEVEL 

POPULATION 

INCIDENCE 

POPULATION 

INCIDENCE 

Low 

10 

3.0 

10 

2.4 

Moderate 

25 

14.3 

70 

60.0 

High 

65 

82.7 

20 

37.6 
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V.  EFFECT  OF  SENSITIVITY  DISTRIBUTION  ON  SKIN  CANCER  INCIDENCE 

It  should  not  be  inferred  from  the  above  example  that  dif¬ 
ferential  sensitivity  is  not  a  very  significant  parameter  in 
skin  cancer  incidence.  To  illustrate  the  large  effect  that  the 
sensitivity  distribution  can  have,  consider  the  case  of  popula¬ 
tions  of  equal  exposure  and  variable  values  of  f^  and  f  ^ .  These 
two  values  define  the  sensitivity  distribution  for  a  population 
divided  into  3  sensitivity  groups  since  fg  ■  1  -  ^  -  f_.  The 
incidence  for  a  given  sensitivity  distribution  is 

1  "  k  [S1  fl  +  S2  +  S3  f3  ]  •  (25) 

where  k  is  a  constant  whose  value  is  dependent  on  exposure. 
However,  it  is  the  relative  incidence  I'  or  I/k  that  is  of  in¬ 
terest  here.  Assuming  the  sensitivity  values  *  1,  “  3, 

and  =  9  of  Case  H, 

I'  =  fl  +  3  f2  +  9  f3  (26) 

*  3  -  2  f1  +  6  f3  .  (27) 

I '  is  plotted  vs  f^  in  Pig.  10,  with  f^  as  a  parameter. 

For  a  given  value  of  f^  the  value  of  f^  is  constrained  by  the 
inequality 

0  s  f  s  1  -  T1  .  (28) 

Note  that  for  a  population  with  a  sensitivity  distribution 
f1  =  0.1,  f^  =  0.6  as  might  be  characteristic  of  a  population 
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FIGURE  10.  Relative  incidence  for  equal  exposure  populations 
vs  sensitivity  distribution. 


dominated  by  residents  of  Scandinavian  or  Irish  descent,  I'  = 
6.2  (point  A  in  Fig.  10);  for  a  population  with  a  distribution 
f1  =  0.3,  f^  =  0.1,  as  might  be  characteristic  of  a  population 
with  many  residents  of  Mediterranean  descent  and  few  of 
Scandinavian  or  Irish,  I'  =  3.0  (point  B  in  Fig.  10).  When  a 
possible  factor  of  2  or  more  in  incidence  due  to  the  effects 
of  differential  sensitivity  is  combined  with  a  possible  factor 
of  2  or  more  due  to  the  effects  of  differential  exposure,  it 
is  not  too  surprising  that  a  factor  of  4  or  more  in  incidence 
could  be  registered  in  different  populations  at  the  same  lati¬ 
tude,  as  illustrated  by  the  data  in  Fig.  1. 
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VI.  EFFECTS  OF  NON-SOLAR- INDUCED  SKIN  CANCER  INCIDENCE 


A  substantial  fraction  of  basal  cell  epitheliomas,  per¬ 
haps  50  percent  in  northern  U.S.  cities  (Ref.  1),  are  not  in¬ 
duced  by  solar  ultraviolet  radiation.  Inclusion  of  such  can¬ 
cers  in  mathematical  models  can  lead  to  significant  errors, 
particularly  when  comparing  incidences  in  different  regions 
to  determine  the  biological  amplification  factor.  This  effect 
is  easily  illustrated  by  making  a  one-dimensional  incidence 
comparison  between  two  regions,  i.e.,  the  sensitivity  of  the 
populations  are  assumed  to  be  equal,  but  aggregate  exposure  is 
represented  by  UV-B  dose  DA  in  region  A  and  Dg  in  region  B. 

If  all  skin  cancers  reported  are  assumed  to  be  induced  by 
solar  ultraviolet  radiation,  the  incidences  in  regions  A  and  B 
would  be  incorrectly  represented  by  the  relations 

IA  -  k  D^0  (29) 

XB  "  k  d/0  ,  (30) 

where  k  is  a  constant  and  nQ  is  the  biological  amplification 
factor,  calculated  on  the  basis  of  the  erroneous  assumption. 
Solving  Eqs .  (29)  and  (30)  for  n  , 


/ 


Since  there  is  an  incidence  I^A  of  non-solar-induced  skin 
cancer  in  region  A  and  an  incidence  Ixg  of  non-solar-induced 
skin  cancer  in  region  B, 


k  da  S+  xxa 

(32) 

n 

k  DB  +  IXB  * 

(33) 

where  n  is  the  amplification  factor  for  solar-induced  skin  can- 
s 

cers  only.  Solving  Eqs.  (32)  and  (33)  for  n  , 

S  i 


n  = 
s 


(34) 


The  error  in  biological  amplification  can  be  expressed  by 
the  ratio  n^nQ  or,  from  Eqs.  (31)  and  (34), 


(35) 


In  Fib.  11,  the  ratio  n  /n  is  plotted  vs  the  non-solar  in- 

s  o 


duced  skin  cancer  incidence  I 
200  and  I 


for  region  B,  for  the  Case  I, 


XB  *  -'-e-'"'  "»  *A 

XA  =  20  or  10  percent  of  IA.  Three  cases  for  Ig  are 

shown  in  Fig.  li:  Ig  =  50,  100,  and  150.  The  values  of  Ixg 
for  each  case  were  assumed  to  lie  between  zero  and  0.5  I.,.  The 

O 


case  n 


nQ  can  only  occur  when  the  fraction  of  non-solar-in- 


duced  skin  cases  is  the  same  in  regions  A  and  B,  or 
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NON-SOLAR-INDUCED  INCIDENCE  IN  REGION  B,  lXB 


FIGURE  11.  Effect  of  non-solar-induced  skin  cancer 

incidence  on  biological  amplification  factor. 
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However,  it  can  be  expected  that  in  general  the  region  of 
significantly  higher  incidence  (region  A)  will  have  associated 
with  it  a  smaller  fraction  of  non-solar  induced  skin  cancers, 
so  that 


and  it  would  follow  that 


s  o 


If  I, 


*  IY.  *  20,  the  error  in  n_  with  respect  to  n  lies 


between  14  percent  (Ig  -  150)  and  30  percent  (Ig  *  50).  If 

IXB  =  1,0»  the  error  in  ns  is  60  percent  for  I  =  100,  and  72 
percent  for  Ig  =  150. 

If  the  Incorrect  biological  amplification  factor  is  used 
to  predict  the  expected  number  N  of  additional  skin  cancer  cases 
that  would  result  in  some  geographic  region  as  a  result  of  a 
hypothetical  increase  in  ultraviolet  dose  AD  following  a  reduc¬ 
tion  in  the  thickness  of  the  stratospheric  ozone  layer,  there 
will  be  an  error  in  N.  If  P  is  the  population  of  the  region, 

the  additional  number  of  skin  cancers  AN  based  on  the  use  of 

o 

biological  amplification  factor  nQ  Is 


-PI 


"o(f)  • 


The  correct  amplification  factor  n  ,  however,  should 

s 

only  be  applied  to  the  solar-related  cases  of  skin  cancer,  or 


4N  -  P  (I  -  Ix>n«(f) 


(40) 


The  error  in  N  may  be  represented  by  the  ratio  AN/AN  or 

«.-(■-  *)  k)  ■ 

In  Pig.  12.  the  ratio  AN/AN,  is  plotted  vs  n  /n  with  the 
non-solar-related  skin  cancer  incidence  Ix  as  a  parameter.  The 
error  in  N  may  be  positive  or  negative.  If  Ix  <  (l-no/nJ  I, 
then  AN  >  ANQ  and  the  expected  number  of  additional  cases  would 
be  underestimated  by  ANq;  if  Ix  >  (1  -  no/ng)  I,  then  AN  <  ANq 
and  the  number  of  additional  cases  would  be  overestimated. 


Ml 
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AN/ANq 


FIGURE  12.  Effect  of  non-sol ar -i nduced  skin  cancers  on  number 
of  predicted  additional  skin  cancers  In  a 
geographic  region. 
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VII.  MODELLING  NON-MELANOMA  SKIN  CANCER  INCIDENCE 
BY  TYPE  OF  TUMOR  AND  SEX 


The  incidence  of  skin  cancer  varies  not  only  with  age,  ex¬ 
posure,  and  sensitivity,  but  also  by  type  of  tumor  and  sex.  The 
total  non-melanoma  skin  cancer  incidence  I  is  the  sex  weighted 
sum  of  six  incidence  terms,  i.e.. 


I  =  M  (I.  +  I  +  I  )  +  (1  -  M)(I.  „  +  I  +  I  ) 

bm  sm  xm  bf  sf  xf 


where  M  is  the  fraction  of  adult  males,  subscripts  m  and  f  denote 
male  and  female,  respectively,  subscripts  b  and  s  denote  basal 
cell  epitheliomas  and  squamous  cell  carcinomas,  respectively,  and 
subscript  x  denotes  non-solar-related  carcinomas.  Each  of  the 
four  solar-related  skin  cancer  incidence  terms  of  Interest  is 
composed  of  the  sum  of  the  incidences  of  each  anatomic  site  which 
vary  in  frequency  with  geographic  region  (Ref.  1).  The  age- 
specific  incidence  of  basal  cell  epithelioma  differs  from  that  of 
squamous  cell  carcinoma  (Ref.  6),  and  the  dose-response  relation¬ 
ship  could  therefore  also  be  expected  to  differ  for  these  two 
types  of  tumors.  Adding  further  to  the  complexity  of  the  skin 
cancer  modelling  problem  is  the  likelihood  that  temperature, 
humidity,  and  wind  speed  affect  skin  cancer  incidence  (Ref.  1). 

It  is  therefore  apparent  that  a  skin  cancer  model,  including  dif¬ 
ferential  exposure  and  differential  sensitivity,  while  much  more 
sophisticated  than  an  aggregate  model,  represents  but  a  first 
step  in  a  comprehensive  quantitative  understanding  of  the 
epidemiology  of  non-melanoma  skin  cancer. 
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VIII.  EFFECTS  OF  CHANGES  IN  STRATOSPHERIC  OZONE 

Most  recent  analyses  of  aircraft  atmospheric  pollution  in¬ 
dicate  that  "the  NOx  constitu“nt  of  engine  exhaust  from  both 
SSTs  at  17  to  20  km  (except  Jn  very  large  numbers)  and  from  sub- 
sonics  at  9  to  14  km,  will  contribute  to  an  increase  in  the 
ozone  column;  water  vapor  emissions,  on  the  other  hand,  at  least 
from  SSTs,  (will)  tend  to  decrease  ozone  unless,  as  some  model¬ 
ing  results  indicate,  radiative  feedback  attempts  compensate. 

The  net  effect  on  the  ozone  column,  of  water  and  NC>x  emissions, 
with  current  N0x  emission  indices,  is  apparently  slightly  posi¬ 
tive  for  both  classes  of  aircraft  in  foreseeable  numbers" 

(Ref.  7). 

Halocarbons  in  the  stratosphere  will  probably  decrease  ozone 
more  than  aircraft  will  increase  ozone.  Halocarbon  effects  and 
aircraft  effects  (and  all  others)  will  change  with  time.  Time- 
dependent  modelling  results  showing  both  effects  are  not  available: 
a  brief  review  of  the  current  literature  suggests  halocarbon  effects 
will  dominate.  However,  stratospheric  ozone  changes  also  may  be 
caused  by  climatic  as  well  as  chemical  changes  (Ref.  8).  For  pur¬ 
poses  of  discussion,  the  effects  of  a  postulated  increase  in  strat¬ 
ospheric  ozone  also  deserve  some  consideration. 

The  magnitude  of  a  decrease  in  solar  UV  radiation  due  to  a 
postulated  increase  in  total  ozone  was  discussed  in  Ref.  9.  The 
relationship  is  symmetrical,  i.e.,  if  a  small  specified  percent¬ 
age  decrease  in  ozone  leads  to  a  factor  K  increase  in  ultra¬ 
violet  exposure  E  or  ultraviolet  dose  D,  an  increase  in  ozone  of 
the  same  percentage  leads  to  a  factor  K-1  decrease  in  D.  For  a 
small  percentage  decrease  in  D,  the  fractional  decrease  in  skin 
cancer  incidence  LI/1  is  given  by  n  AD/D.  The  skin  cancer 
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mathematical  model  itself  is,  of  course,  unaffected  by  the  di¬ 
rection  of  the  change  in  the  input  parameter.  The  biological 
consequences  of  a  thicker  ozone  layer  are  another  matter.  Al¬ 
though  skin  cancer  incidence  would  decrease  —  an  obvious  benefit 
to  society  —  the  healthful  effects  of  solar  ultraviolet  radia¬ 
tion  would  also  be  reduced.  These  healthful  effects  include 
prevention  of  rickets  and  beneficial  skin  effects  which  do  not 
readily  lend  themselves  to  quantification.  Among  these  are 
skin  conditions  such  as  psoriasis,  sebhorric  dermatitis,  acne, 
and  actopic  exema,  conditions  in  which  all  practicing  dermatolo¬ 
gists  note  an  improvement  of  the  skin  during  the  summer  months 
when  exposure  to  solar  radiation  is  increased.* 


Personal  communication  with  Dr.  T.  P.  Nigra,  Chairman, 
Dermatology  Section,  Washington  Hospital  Center. 
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IX.  IMPLICATIONS  OF  MODEL  RESULTS  TO  PREDICTION  OF 
INCREASE  IN  SKIN  CANCER  INCIDENCE 


The  technique  most  widely  used  to  predict  the  increase  in 
skin  cancer  incidence  resulting  from  a  hypothetical  decrease 
in  the  thickness  of  the  stratospheric  ozone  layer  is  to  compare 
incidence  values  for  geographic  regions  of  different  latitude 
with  calculated  values  of  the  corresponding  UV-B  dose  incident 
on  the  ground  (e.g.,  Ref.  10).  Arbitrary  adjustment  is  made  for 
aggregate  population  exposure  differences  due  to  differences  in 
length  of  season  and  climate,  an  adjustment  which  was  unneces¬ 
sary  in  the  comparisons  investigated  in  this  paper  because  New 
York  City  and  Iowa  are  at  approximately  the  same  latitude.  The 
results  of  the  model  illustrate  that  a  refinement  of  this  tech¬ 
nique  should  include  the  important  effects  of  differential  ex¬ 
posure  and  differential  sensitivity  within  the  population  of 
each  of  the  geographic  regions  used  for  comparison  purposes. 
Clearly,  such  a  refinement  is  difficult  if  not  impossible  to 
implement  with  the  present  lack  of  knowledge  of  the  distribution 
of  population  exposure  and  the  distribution  of  sensitivity  in 
any  given  geographic  region.  It  is  to  be  hoped  that,  sometime 
in  the  near  future,  there  will  exist  a  data  bank  of  population 
exposure  distributions  based  on  surveys  with  personal  DUV  dosim¬ 
eter  readings  (such  dosimeters  are  currently  in  research  and 
development*).  There  already  exists  a  DUV  network  which  is  pro¬ 
viding  information  on  the  dally  DUV  dose  incident  on  the  ground 


One  such  personal  dosimeter  is  currently  under  development 
by  Ricker  Laboratories  Inc.,  St.  Paul,  Minnesota. 
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for  each  day  of  the  year  (Ref.  11).  There  are  also  dermatolo¬ 
gists  investigating  the  difficult  problem  of  finding  a  suitable 
measure  of  individual  sensitivity  to  skin  cancer  based  on  char¬ 
acteristics  of  the  skin,  e.g.,  time  to  show  a  minimal  erythemal 
dose*.  A  complicating  factor  is  the  effect  of  tanning  and  pig¬ 
mentation.  Long-term  exposure  to  sunlight  produces  tanning  in 
some  individuals  and  thereby  reduces  their  sensitivity  to  light. 
An  additional  approach  could  be  genetic  predisposition  to  skin 
cancer.  Techniques  are  at  hand  to  genetically  type  individuals 
with  regard  to  their  HLA,  B,  C,  and  D  loci  (tissue  typing  used 
in  transplant  immunology)**.  If  a  sensitivity  measure  can  be 
found  and,  with  more  detailed  and  accurate  epidemiological  data 
than  are  now  available,  demonstrated  to  be  satisfactorily  corre¬ 
lated  with  skin  cancer  incidence,  the  necessary  foundations  will 
have  been  laid  for  the  inclusion  of  differential  exposure  and 
differential  sensitivity  into  a  model  predicting  increases  in 
skin  cancer  incidence  that  would  follow  a  stratospheric  ozone 
reduction.  Such  a  model  should  incorporate  the  methodology  out¬ 
lined  in  this  paper,  but  the  number  of  exposure-sensitivity 
groups  should  probably  be  Increased  from  9  to  16  or  more. 


Such  as  the  Wucherpf ennig  method  of  measuring  minimum 
erythemal  dose  which  consists  of  spatially  filtered  wheels 
providing  varying  levels  of  exposure  (described  at  a  meet¬ 
ing  sponsored  by  the  Association  Internationale  de  Photo- 
biologie  in  Lausanne,  Switzerland,  September  26-29,  1978). 


«* 


Personal  Communication  with  Dr.  T.  P.  Nigra,  Chairman, 
Dermatology  Section,  Washington  Hospital  Center. 


X.  CONCLUSIONS 


The  principal  conclusions  that  can  be  drawn  from  the  anal¬ 
ysis  in  this  paper  are  the  following: 

1.  Differential  exposure  and  differential  sensitivity  in 
a  population  have  a  significant  impact  on  skin  cancer 
Incidence,  and  should  be  taken  into  consideration  when 
comparing  skin  cancer  incidence  in  populations  resid¬ 
ing  in  different  geographic  regions. 

2.  It  is  possible  to  mathematically  model  skin  cancer 
incidence,  including  the  effects  of  differential 
exposure  and  differential  sensitivity. 

3.  Parameters  identified  as  essential  to  the  formulation 
of  a  mathematical  model  are  the  following: 

a.  The  range  of  exposure  values  E1> 

b.  The  probability  density  function  of  exposure 
P  (E1). 

c.  The  range  of  sensitivity  values  Si> 

d.  The  probability  density  function  of  sensitivity 
p  (S^  or 

e.  The  biological  amplification  factor  n. 

4.  A  mathematical  model  of  skin  cancer  incidence  based  on 
a  division  of  the  population  into  three  exposure  and 
three  sensitivity  groups  was  derived  and  applied  to  the 
case  of  New  York  City.  A  set  of  22  cases  was  investigated. 
Exposure  and  sensitivity  values  which  appear  to  best  fit 
the  input  incidence  data  were  determined  by  eliminating 


those  cases  which  led  to  negative  probabilities.  The 
remaining  solutions  were  then  subjected  to  criteria 
designed  to  minimize  departures  from  expected  behavior. 
Of  the  cases  investigated,  the  best  New  York  City  solu¬ 
tion  had  the  following  set  of  values:  n  =  1,  E^  =  1, 

E2  =  4,  E3  -  8,  Sx  •  1,  S2  •  3,  S3  =  9,  and  r1  -  f = 

0.1. 

5.  A  reversal  of  the  New  York  City  population  exposure 
distribution  was  assumed  for  a  rural  region  such  as 
Iowa,  assuming  no  change  in  the  population  sensitivity 
distribution.  It  was  found  that  Iowa  could  have  more 
than  double  the  incidence  of  New  York  City  under  the 
assumed  conditions. 

6.  Two  white  populations  with  markedly  different  sensi¬ 
tivity  distributions  and  equal  exposure  could  have 
incidences  that  differ  by  a  factor  of  2  or  more. 

7.  The  effect  of  not  eliminating  non-solar-related  skin 
cancers  in  skin  cancer  incidence  data  on  the  prediction 
of  increased  incidence  due  to  stratospheric  ozone  re¬ 
duction  is  to  introduce  an  error  which  may  be  positive 
or  negative. 

8.  In  order  to  effectively  utilize  the  mathematical  model 
derived  in  this  paper  for  predicting  skin  cancer  in¬ 
creases  resulting  from  stratospheric  ozone  depletion, 
it  will  be  necessary  to  pursue  the  following  programs: 

a.  Develop  a  personal  UV-B  dosimeter  (in  progress). 

b.  Conduct  exposure  surveys  of  individuals  wearing 
a  personal  UV-B  dosimeter  in  various  population 
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groups  residing  in  various  geographic  regions* 
for  which  accurate  skin  cancer  incidence  data 
is  available. 


c.  Determine  a  measure  of  sensitivity  to  skin  can¬ 
cer  which  can  be  correlated  to  incidence  and 
which  is  easily  applied  in  the  field  (under 
investigation*  * ) . 

d.  Classify  various  population  groups  residing 
in  various  geographic  regions  for  which  skin 
cancer  incidence  data  is  available  according 
to  (c). 

e.  Pathology  criteria  are  needed  to  distinguish 
skin  cancers  that  are  caused  by  solar  ultra¬ 
violet  radiation  from  those  that  are  not. 

f.  Non-melanoma  skin  cancer  incidence  is  cur¬ 
rently  reported  to  cancer  registries  by  age 
and  sex.  It  should  also  be  reported  by  type 
of  tumor  (basal  cell  epithelioma  and  squamous 
cell  carcinoma),  by  relation  to  solar  ultra¬ 
violet  radiation,  if  and  when  item  (e)  above 
becomes  feasible,  and  by  sensitivity,  if  and 
when  item  (c)  becomes  feasible. 


In  England,  for  example.  Dr.  I.  Magnus  has  reported  (Lausanne, 
Switzerland,  Sept.  1978)  that  tests  with  his  personal  UV 
badges  indicate  that  sendentary  people  receive  approximately 
three  percent  of  the  ambient  DUV  radiation. 


Dr.  J.  C.  Van  der  Leun  has  suggested  (Lausanne,  Switzerland, 
Sept.  1978)  measures  of  carcinogenesis  susceptibility  other 
than  minimum  erythemal  dose  such  as  increase  in  redness  with 
increase  in  dose,  and  an  erythema  index  which  measures  the 
increase  in  the  size  of  the  red  area  with  an  increase  In 
dose. 
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